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Abstract 

Noninvasive  functional  imaging,  monitoring  and  quantification  of  microbubbles  forming  in 
blood  and  tissues  upon  rapid  changes  in  barometric  pressure  are  extremely  important  for 
effective  therapy  and  diagnostics  of  several  diseases  as  well  as  for  imaging  and  drug  delivery 
projects.  However,  current  techniques  are  unable  of  imaging  and  efficient  detection  of  bubbles 
with  diameter  less  than  50  micrometers.  The  goal  of  this  proposal  was  to  develop  novel  Phase- 
Sensitive  Swept  Source  Optical  Coherence  Tomography  (PhS-SSOCT)  technique  capable  of 
real-time,  sensitive,  accurate,  and  noninvasive  imaging,  monitoring,  and  quantification  of 
microbubbles  in  tissues.  During  these  studies,  a  novel  phase  resolved  system  based  on  Swept 
Source  Optical  Coherence  Tomography  (SSOCT)  has  been  developed.  The  system  has  an  axial 
resolution  of  10  pm,  phase  sensitivity  of  0.03  radians,  imaging  depth  of  up  to  6  mm  in  air,  and 
in-depth  scanning  speed  of  20  kHz  for  a  single  A-line.  The  performance  of  the  sensing  system 
was  carefully  evaluated  in  optical  phantoms  containing  gas  microbubbles  with  different 
diameter.  Obtained  results  demonstrate  that  bubbles  with  diameter  greater  than  10  pm  could  be 
detected  by  both  structural  imaging  and  phase  response  whereas  bubbles  with  diameters  less  than 
10  pm  could  be  detected  by  the  phase  response  of  the  SSOCT  with  high  sensitivity.  The 
accuracy  for  measurement  of  the  diameter  of  gas  microbubbles  is  limited  to  10  pm  in  structural 
imaging  and  0.01  pm  in  phase-sensitive  monitoring.  Preliminary  studies  were  also  performed  in 
animals  in  vivo  for  the  rapid  assessment  of  the  circulating  microbubbles.  The  in  vivo  results 
demonstrate  capability  of  the  developed  instrument  to  image  bubbles  with  diameter  above  8  pm. 
The  results  from  the  study  indicate  that  the  developed  SSOCT  system  could  be  used  to  detect 
fast-moving  microbubbles  and  warrant  further  investigation  in  biological  tissues  in  vivo. 
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Introduction 


Formation  of  microbubbles  in  the  human  body  has  been  extensively  studied  by  many 
researches  since  discovery  of  Caisson  disease  by  Paul  Bert  in  1878.  In  its  classic  and  most  severe 
form,  gas  emboli  present  catastrophically  (~4%  of  victims)  with  collapse,  loss  of  consciousness, 
apnea,  and  cardiac  arrest.  Formation  of  microbubbles  inside  the  body  could  be  attributed  due  to 
different  physiological,  physical,  and  intervention  mechanisms.  These  bubbles  can  travel  to  any 
part  of  the  body,  accounting  for  many  serious  (and  sometimes  life-threatening)  disorders  [1],  Gas 
bubbles  in  the  back  or  joints  can  cause  localized  pain  (the  bends).  In  the  spinal  cord  or  peripheral 
nerve  tissues,  bubbles  may  cause  paresthesias,  neurapraxia,  or  paralysis.  A  bubble  forming  in  the 
circulatory  system  can  lead  to  pulmonary  or  cerebral  gas  emboli.  Formation  and/or  introduction 
of  gas  microbubbles  in  human  blood  and  tissues  remains  a  serious  long-term  sequel  in  patients 
undergoing  cardiac  valve  replacement  (with  an  annual  risk  of  up  to  4%)  [2-4];  cardiopulmonary 
bypass  and  other  open-heart  surgeries  [3-8];  high-intensity  focused  US  therapy  [9-11];  cesarean 
section  and  operative  hysteroscopy  [12-24];  orthopedic  surgery  [25,  26];  and  various  laser 
ablation  and  laparoscopic  surgeries  [27-30],  Additionally,  gas  embolism  happens  in  endoscopy 
[31],  tissue  biopsy  [32],  neurosurgery  [33],  liver  transplantation  [34-36],  during  central  venous 
line  insertion  and  removal  [37,  38],  and  even  during  intravenous  antibiotic  delivery  at  home  [39], 
The  use  of  ultrasound  bubble  contrast  media  could  also  lead  to  the  emboli  [40,  41].  Moreover, 
FDA  and  Bristol  Myers  Squibb  Imaging  issued  a  new  alert  and  important  safety  changes  to  the 
prescribing  information  for  DEFINITY®  (Perflutren  Lipid  Microsphere,  a  diagnostic  drug  used 
as  a  contrast  enhancement  during  echocardiographic  procedures)  in  October  2007  [40],  A  new 
clinical  study  found  serious  cardiopulmonary  problems  including  fatalities  associated  with  the 
use  of  DEFINITY®  microbubbles  likely  attributed  to  venous  gas  emboli.  The  maximum  diameter 
of  the  DEFINITY®  microbubbles  is  20  pm.  Thus,  formation  and  introduction  of  gas 
microbubbles  in  human  blood  and  tissues  is  not  only  a  problem  for  the  divers  but  also  is  a 
significant  everyday  clinical  problem  affecting  thousands  of  patients  undergoing  various  surgical 
and  therapeutic  procedures. 

Bubbles  may  act  as  emboli  and  block  circulation,  as  well  as  cause  mechanical  compression 
and  stretching  of  the  blood  vessels  and  nerves  [42-45].  Additionally,  the  blood-bubble  interface 
acts  as  a  foreign  surface,  activating  the  early  phases  of  blood  coagulation  and  the  release  of 
vasoactive  substances  from  the  cells  lining  the  blood  vessels  and  may  trigger  activation  of 
inflammatory  response  and,  therefore,  expression  of  different  proteins  and  protein  patters  [46- 
54],  This  can  further  worsen  gas  emboli  symptoms.  However,  the  underlying  pathophysiology  is 
still  poorly  understood  [55],  Currently,  there  is  no  reliable  method  for  prediction  or  diagnosis  of 
bubble-associated  diseases  prior  to  commencement  of  apparent  clinical  symptoms  (e.g.  severe 
pain  in  joints,  pulmonary  problems,  disorientation  and  mental  dullness,  vomiting,  skin  rush,  and 
coma).  It  is  widely  recognized  that  prompt  detection  of  microbubbles  in  tissues  and  blood  is 
the  key  for  successful  management  and  treatment  of  persons  with  different  forms  of  gas 
emboli  [1].  Early  noninvasive  detection  of  bubbles  forming  in  tissues  is  required  for  effective 
treatment  and  prediction  of  severity  of  emboli  as  well  as  for  many  therapeutic  and  research 
applications  utilizing  bubbles  as  contrast  or  drug  delivery  agents. 

Depending  on  the  clinical  situation,  the  nature  of  the  gas  emboli  and  the  number  of  embolic 
events  can  vary  greatly.  For  example,  Hills  and  Butler  [56]  measured  intravascular  gaseous 
emboli  ranging  from  19  up  to  700  pm  following  decompression  in  living  dogs.  Gersh  detected 


bubble  sizes  between  60  and  300  jam  in  both  intravascular  and  extravascular  sites.  However, 
several  studies  suggested  that  bubbles  with  diameter  as  small  as  8  pm  could  cause  blockage  and 
result  in  the  trauma  and  onset  of  the  symptoms  [40,  57,  58].  Therefore,  in  order  to  be  effective, 
an  imaging  or  sensing  technique  should  accurately  detect  bubbles  with  diameter  >8  pm. 

The  number  of  embolic  signals  can  also  show  wide  variations.  Georgiadis  et  al.  recorded 
between  0  and  620  embolic  events  per  30  min  period  for  patients  with  prosthetic  heart  valves 
[59],  Mullges  et  al.  have  observed  between  0.53  and  59.05  embolic  signals  per  minute  during 
extracorporeal  circulation  in  patients  undergoing  cardiac  surgery  [60].  However,  these 
observations  were  limited  by  the  resolution  of  utilized  technique  (that  is  50  pm  at  the  best)  and 
actual  number  of  “trouble”  microbubbles  (diameter  >  8  pm)  remains  unknown.  Accordingly, 
development  of  a  device  that  allows  detection  of  small  microbubbles  with  sufficient  accuracy 
and  sensitivity  would  be  a  valuable  tool  for  embolism  studies,  either  for  clinical  or  pre-clinical  or 
both. 

Previously,  several  imaging  techniques  have  been  proposed  and  applied  to  study 
microbubbles  in  blood  and  tissues  including  Doppler  Sonography,  MRI,  Nuclear  Imaging,  and 
Computer  Tomography.  Doppler  Sonography,  the  most  popular  technique  since  air-bubble 
interface  produce  strong  ultrasonic  reflection  in  the  region  of  1  MHz  to  20MHz,  is  an  ultrasound 
diagnostic  imaging  technique,  enhanced  with  Doppler  effect,  capable  of  assessment  moving 
bubbles  by  calculating  frequency  shift  of  a  particular  sample  volume  [61,  62],  However,  Doppler 
Sonography  can  detect  only  moving  intravascular  bubbles  with  a  diameter  of  approximately  50 
pm  [63-65].  Evidently,  the  resolution  of  this  imaging  method  should  be  improved  in  order  to 
achieve  sensitive  imaging  and  assessment  of  small  pm-sized  bubbles.  Also,  new  method  has  to 
be  developed  for  detection  of  stationary  bubbles  in  tissue  that  could  be  used  for  prediction  or 
diagnosis  of  gas  emboli  [66], 

Optics-based  techniques  have  great  intrinsic  potential  to  achieve  the  goal  of  noninvasive 
imaging  of  microbubbles  in  epithelial  tissues  and  blood  microvessels.  Confocal  laser  scanning 
microscopy  (CLSM)  [67],  two-photon  fluorescence  microscopy  (2P-FM)  [68-72],  and  higher 
harmonic  generation  (HHG)  microscopy  [73,  74]  are  some  examples  of  optical  methods  that 
have  been  applied  in  different  fields  of  biological  research.  CLSM  have  significant  axial  and 
lateral  resolutions,  but  it  is  limited  to  ~100  pm  penetration  depth  due  to  high  attenuation  of 
visible/ultraviolet  excitation  light.  2P-FM  has  higher  penetration  depth  due  to  near-infrared 
(NIR)  excitation  of  fluorophores,  but  general  needs  of  exogenous  fluorophores  make  this 
technique  not  truly  “noninvasive”.  HHG  microscopy  does  not  require  application  of  exogenous 
fluorophores,  but  is  very  costly  (a  commercial  Ti:sapphire  system  with  just  the  basics  goes  for 
over  $150K)  and  bulky  [75], 

Optical  interferometric  techniques  are  extremely  sensitive  to  local  changes  in  scattering, 
absorption,  and  refractive  index  of  the  tissues  and  cells.  Since  the  average  refractive  indexes  of 
blood,  skin,  and  air  are  quite  different  (1.4,  1.55  and  1.0,  respectively  in  NIR),  an  optical-based 
sensor  will  be  capable  of  assessing  formation  of  gas  bubbles  with  ultra-high  sensitivity  and 
accuracy  [76,  77].  Here  we  describe  our  progress  in  development  of  a  novel  functional  biosensor, 
based  on  phase-sensitive  Optical  Coherence  Tomography  (OCT)  technique,  for  noninvasive, 
accurate,  and  sensitive  imaging,  monitoring,  and  assessment  of  microbubbles  in  skin  and  skin’s 
blood  mirovessels  by  using  both  amplitude  imaging  and  phase-resolved  sensing  methods. 


OCT  is  a  relatively  novel  non-invasive  optical  diagnostic  technique  that  provides  depth- 
resolved  images  of  tissues  with  resolution  up  to  few  micrometers  at  depths  of  up  to  several 
millimeters.  This  technique  was  introduced  in  1991  to  perform  tomographic  imaging  of  the 
human  eye  [78],  Since  then,  OCT  is  being  actively  developed  by  several  research  groups  mainly 
for  clinical  diagnostic  applications  (reviewed  in  [79-81]). 

Current  OCT-based  imaging  techniques  may 
be  divided  into  two  classes:  time  domain  OCT 
(TDOCT)  and  spectral  domain  OCT  (SOCT)  (Figure 
1).  While  conventional  TDOCT  methods  are  widely 
used  in  clinical  and  research  laboratories,  recently,  the 
significant  signal-to-noise  (SNR)  advantage  of  SOCT 
over  TDOCT  (up  to  30  dB)  has  been  demonstrated 
[82-84],  Also,  simplicity  of  the  interferometer 
construction,  the  absence  of  mechanical  in-depth 
scanning,  and  the  SNR  advantage  of  SDOCT  are 
facilitating  its  further  development  for  many 
applications.  Moreover,  the  relaxed  requirement  for 
mechanical  in-depth  scanning  in  SOCT  systems 
allows  using  a  common  beam  path  for  the  reference 
and  sample  arms  of  the  interferometer  with  the 
reference  plane  located  near  the  surface  of  the  sample 
[85].  Such  a  common-path  interferometer  is  easier  to 
align  and  is  less  sensitive  to  mechanical  vibrations, 
and  it  allows  for  more- straightforward  compensation 
for  dispersion  and  polarization  differences  between 
reference  and  sample  arms. 

SOCT  may  be  further  divided  into  two  distinct 
methods  of  acquiring  spectral  information  from  an 
object  under  study:  1)  Fourier-domain  OCT  (FDOCT) 
-  using  a  broadband  laser  source  (typical  FWHM  >75 
nm)  with  a  grating  and  a  photodiode  array  in  the 
detection  arm;  and  2)  swept-source  OCT  (SSOCT)  -  using  a  rapidly  tunable,  narrow-linewidth 
laser  source  over  a  broad  optical  bandwidth  with  a  single  photodiode  in  the  detection  arm.  While 
FDOCT  and  SSOCT  have  similar  sensitivity  advantages  over  conventional  TDOCT  systems 
[82],  the  SSOCT  system  can  detect  interference  fringes  over  a  substantially  longer  range  of  time 
delays  between  reflections  from  reference  and  sample  interfaces. 

Phase-sensitive  measurements  of  biological  objects  provide  additional  information  on  optical 
properties  and  play  a  key  role  in  development  of  functional  OCT  (fOCT)  methods.  Recently 
introduced  TDOCT  with  differential  phase  contrast  [86,  87]  shows  great  potential  for  imaging 
and  classification  of  tissues  as  well  as  for  fOCT  [88].  While  conventional  OCT  is  based  on  the 
detection  and  analysis  of  amplitude  (or  spectrum)  of  back-scattered  optical  radiation,  the  phase 
sensitive  OCT  utilizes  the  phase  information  obtained  from  two  tissue  interfaces  by  the  same  or 
two  probing  beams.  Variations  in  the  sample  refractive  index  will  be  reflected  in  the  phase 

difference,^,  between  these  two  signals  from  different  interfaces.  The  phase  sensitive  OCT 
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Figure  1:  Schematic  diagram  of  time  domain 
(upper  part)  and  spectral  domain  (lower  part) 
OCT  systems. 


technique  (that  could  be  both  TDOCT  or  FDOCT)  is  capable  of  measuring  Angstrom  scale  path 

Ap  =  —A(p 

length  change  between  the  beams  (associated  with  the  phase  difference  as  4n  )  in  tissues 
[88,89], 

In  OCT,  a  two-beam  interferometer  is  used  to  detect  the  backscattered  photons  from  a  tissue 
of  interest  within  the  coherence  length  of  the  laser  source.  Typically,  a  broadband  low  coherence 
laser  source  is  split  into  two  arms  (reference  arm  and  sample  arm).  The  scattered  light  from  the 
sample  is  captured  in  such  a  way  that  it  recombines  with  the  light  from  the  reference  arm  which 
results  in  the  formation  of  fringes,  provided  the  distance  travelled  by  light  in  both  arms  does  not 
differ  by  more  than  the  coherence  length  of  the  laser  source.  The  interference  signals  are 
recorded  by  the  photo  detector  which  is  given  by  the  Equation  1 : 

i(t)  oc  2  JP^PS  cos  2  k8(t)  ,  (1) 

where  i(t)  is  the  recorded  signal  from  the  detector  and  S(t)  is  the  introduced  path  difference  . 
Pr,  Ps  are  the  optical  powers  from  the  reflected  and  the  sample  arms,  respectively,  and  k  is  the 
wave  number.  The  equation  shows  only  the  signal  in  the  interferometric  term  neglecting  zero 
frequency  component.  Rewriting  equation  1  will  give: 

i(t)  a  Real  {/  5(a>  -  ,  (2) 

where  A0(oo)  =  ^  *  Ad  (&>)  is  the  phase  delay  due  to  the  path  difference  S  (co),  and  S(a>)  is  the 
spectrum  of  the  source.  The  i(t)  and  Aafco)  are  Fourier  Transform  pairs,  thus  by  applying  a 

2n 

complex  Fast  Fourier  Transform  (FFT)  algorithm  to  i(k),  which  is  i(t)  in  k-space  ( k  =  — ),  one 

can  obtain  a  1-D  OCT  depth  profde.  By  generating  transversal  set  of  similar  1-D  depth  profiles 
with  a  scanning  galvo-mounted  mirror,  a  2-D  image  could  be  constructed.  Thus,  2D  image 
contains  both  axial  and  transverse  information. 

For  a  Swept  Source  laser,  k  (t)  =  2n/X(t)  is  the  wave  number  that  does  not  always  obey  a 
linear  relation  k(t)  =  ko+  kit.  In  theory  if  it  obeys  a  linear  relation,  then  the  laser  would  have 
effectively  mapped  the  k-space  to  time  domain,  i.e.,  i(t)  =  i(k).  In  reality,  the  frequency  sweep  of 
the  chirped  laser  is  non-linear.  The  k(t)  contains  higher  order  terms  and  causes  non-linearity  in 
the  frequency  sweep  which  leads  to  a  non-uniform  sampling  interval.  As  the  Fourier  transform 
is  applied  to  get  the  depth  information,  uniformly-spaced  samples  are  required.  By  performing 
Non-Uniform  Fourier  Transformation  (NUFT),  depth  information  from  the  non-uniformly 
spaced  samples  can  be  obtained.  Another  approach  is  by  using  the  interference  fringes  from  a 
fabry-perot  interferometer,  uniformly  spaced  samples  are  obtained.  The  non-linearity  in 
frequency  sweep  is  analyzed  by  several  groups  and  proposed  different  methods  to  overcome  the 
problem.  In  this  design  we  employs  a  Mach-Zehnder  interferometer  based  optical  clock  (MZI- 
OC)  which  generates  equally-spaced  frequency  interferogram  that  range  from  12.5  GHz  to  200 
GHz.  All  peaks,  as  well  as  the  zero  crossings  in  the  recorded  fringes  of  MZI-OC,  are  always 
equally  spaced  in  optical  frequency  space  and  are  used  for  recalibration.  The  samples  are 
collected  at  these  zero  crossings  thus  eliminating  the  non-linearity  in  the  sampling  interval.  This 
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Figure  2:  Improvement  in  resolution  of  SSOCT  after 
calibration 


also  serves  the  purpose  of  converting  into 
the  uniform  k-space.  Without  calibrating  and 
remapping  to  uniform  k-space,  resolution 
degradation  is  observed  along  the  z-space.  In 
addition  to  that,  this  correction  offered  a 
significant  improvement  in  the  resolution  at  a 
particular  depth  as  shown  in  Figure  2:  the 
reflectivity  profile  from  a  weak  reflector 
without  calibrating  (dotted)  and  with  calibrating 
(solid)  of  the  signal. 

In  this  report,  we  present  results  on 
development  of  a  phase  resolved  sensing  system 
based  on  SSOCT  and  evaluated  its  performance 
by  quantifying  microbubbles  of  different 
diameters  in  clear  and  scattering  tissue- 


simulating  media  and  in  pilot  in  vivo  experiments. 


Development  of  Phase-Sensitive  Swept  Source  Optical  Coherence  Tomography  System 


SLS 

(1315+j,-55nrn> 


A  schematic  diagram  of  the  developed  Phase-Sensitive  SSOCT  (PhS-SSOCT)  system  is 
shown  in  Figure  3.  The  system  consists  of  four  main  units:  the  source,  the  interferometer  (Mach- 

Zehnder  in  this  particular 
setup),  the  calibration 
system  and  the  data 
acquisition  electronics.  The 
laser  source  output  is  split 
into  two  arms:  one  arm 
containing  the 

interferometer  and  the  other 
arm  containing  the 
recalibration  and  triggering 
unit.  A  90-10  fiber  coupler 
(Thorlabs)  is  used  to  send 
90%  of  the  light  to  the 
interferometer  so  that  the 
sample  arm  gets  the 
maximum  power.  The 
remaining  10%  is  further 
split  by  a  99-1  fiber  coupler 
with  99%  going  to  a  Fiber 
Bragg  Grating  (FBG)  and 
1%  to  the  MZI-OC.  The 
99%  light  goes  to  the  FBG 
through  a  three-arm 
circulator  and  the  reflected 


Figure  3:  Basic  set  up  of  Ph-SSOCT  for  imaging  of  microbubbles  (ADC: 
Analog  to  Digital  Converter,  BPD:  Balanced  Photo  Detector,  C:  Collimators, 
FBG:  Fiber  Bragg  Grating,  G:  Galvanometer  Scanner,  L:  Lens,  RM: 
Reflective  Mirrors,  SLS:  Swept  Laser  Source,  TTL:  TTL  Signal  Generator). 


pulse  is  passed  to  the  detector  by  the  third  arm  of  the  circulator.  This  detector  outputs  a  voltage 
signal  which  is  converted  to  an  electrically  tunable  TTL  pulse  by  a  pulse  generator  (Stanford 
Research  Systems,  Inc.).  The  TTL  signal  is  tuned  to  a  required  duty  cycle  and  used  to  trigger  the 
analog  to  digital  converter  (ADC).  The  other  1%  is  fed  to  the  MZI-OC  whose  signal  is  detected 
in  balanced  detection  mode,  and  these  electric  signals  from  the  detector  are  acquired  by  one  of 
the  channels  of  the  ADC.  In  the  interferometer  arm,  the  90%  light  is  split  into  1%  and  99%,  each 
going  to  the  reference  arm  and  sample  arm  respectively  via  circulators.  The  light  coming  from 
the  reference  arm  is  passed  through  an  adjustable  pin  hole  to  allow  attenuation  if  required.  The 
reflected  light  from  the  reference  arm  and  sample  arm  are  coupled  into  a  50-50  fiber  coupler 
where  they  recombine  and  form  the  interference  fringes.  These  fringes  are  then  detected  by  a 
balanced  photo  detector  (BPD),  which  subtracts  the  two  signals  to  remove  the  common  mode 
noise.  As  the  fringes  would  be  out  of  phase,  this  BPD  is  effectively  adding  the  fringes  but 
subtracting  the  common  mode  noise.  After  removing  the  common  mode  noise,  the  fringe 
encoded  voltage  is  then  amplified  by  a  Transimpedance  amplifier  (TIA)  and  then  RF  modulated 
and  acquired  by  a  PC  through  the  other  channel  of  the  ADC.  Both  the  information  from  MZI  and 
the  interferometer  is  acquired  simultaneously  by  the  ADC  with  the  receiving  of  the  trigger  from 
the  FBG. 

The  fringes  are  acquired  by  a  14-bit  high-speed  digitizer  (PCI  5122,  National  Instruments). 
The  digitizer  is  operated  at  50  MS/s  and  acquires  2500  sample  points  per  A-scan.  The  sampling 
rate  is  chosen  by  considering  the  fact  that  the  sampling  interval  should  be  smaller  than  the 
instantaneous  line-width  (0.1  nm);  otherwise,  a  large  sensitivity  drop  off  along  the  depth  scan 
would  be  observed.  Out  of  these  2500  points,  the  first  200  and  last  252  points  are  deleted  so  as 
to  select  the  data  corresponding  to  the  laser  wavelength  scans.  Both  the  raw  signal  and  the  MZI- 
OC  signal  will  now  have  2048  points.  After  recalibration  using  MZI-OC,  the  number  of  points  is 
decreased  due  to  the  fact  that  the  number  of  peaks  and  zeros  are  always  less  than  the  total 
number  of  points.  As  a  matter  of  fact,  if  the  nearest  neighborhood  algorithm  is  used  to  find  the 
peaks,  which  gives  1  point  out  of  3  points,  only  one-third  of  the  sample  signal  would  be  utilized. 
Thus,  the  number  of  peaks  and  zeros  registered  are  around  600  indicating  the  presence  of  600 
raw  data  points.  Using  spline  interpolation,  three  points  are  inserted  between  consecutive  raw 
data  points  thus  giving  2400  points.  Again,  the  signal  is  windowed  to  have  2048  points.  It  must 
be  noted  that  the  points  that  are  deleted  from  the  raw  signal  should  not  contain  any  information 
from  the  signal.  If,  in  the  process  of  getting  2048  from  2400,  any  information  containing  the 
fringes  is  lost,  then  instead  of  getting  2048  points  in  the  first  deletion,  more  points  are  retained. 
The  selection  of  points  is  made  in  such  a  way  that  the  calibrated  signal  always  has  2048  points 
without  losing  any  fringe  information.  A  complex  FFT  of  this  signal  gives  the  depth  profile  and, 
due  to  the  symmetry  of  FFT,  each  A-line  corresponding  depth  is  constructed  using  1024  points. 
It  is  worth  noting  that  inserting  a  lower  number  of  points,  e.g.  two,  might  result  in  losing  some  of 
the  higher  frequency  components  which  in  turn  might  reduce  the  imaging  depth  of  the  system. 
Since  the  back-scattered  light  at  a  depth  greater  than  the  coherence  length  of  the  laser  source 
cannot  form  fringes,  the  maximum  imaging  depth  would  not  exceed  the  coherence  length  of  the 
laser  source. 

Signal  processing  includes  several  steps,  including  reference  subtraction,  recalibrating  and 
resampling  to  uniform  k-space  and  image  construction.  The  background  signal  is  recorded  at  the 
beginning  of  every  acquisition  by  blocking  the  sample  arm.  This  signal  contains  the  1%  residual 
light  reflected  from  the  reference  arm  and  from  the  residual  signal  in  the  detectors  due  to 


imperfect  symmetry  of  the  BPD.  Subtracting  the  reference  from  the  interference  signal  helps  to 
increase  the  contrast  and  remove  the  artifacts.  It  removes  the  low-frequency  components 
introduced  by  the  reference  or  ambient  light.  Each  depth  profile  is  obtained  from  the  Fourier 
transform  of  the  interference  fringe  signal.  The  logarithm  of  the  absolute  value  of  this  complex 
FFT  is  mapped  to  the  gray  scale  to  get  an  image  of  a  single  A-scan.  By  generating  a  transversal 
set  of  similar  1-D  depth  profiles  with  the  scanning  galvo-mounted  mirror,  a  2-D  image  is 
constructed.  Thus,  the  2-D  image  contains  both  axial  and  transverse  information.  By  scanning 
the  Y-mirror,  enface  imaging  is  achieved,  thus  allowing  for  3-D  imaging. 
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Figure  4:  (a)  FBG  analog  output  and  the  raw  signal  (b)  oscilloscope  trace  showing  the  synchronization  of  FBG 
trigger,  raw  signal  and  MZI-OC  (c)  the  temporal  phase  response  of  500  pm  cuvette  at  500  pm  depth  before 
stabilization  and  (d)  after  stabilization  (note  the  scale) 

Jitter  present  on  the  output  source  spectrum  causes  the  acquired  fringe  data  to  slowly  drift 
over  time.  These  drifts  introduce  varying  delays  between  the  trigger  signal  and  the  subsequent 
digital  fringe  data.  As  the  phase  is  extracted  from  the  complex  FFT  of  the  fringe  signal,  the 
mismatch  in  fringe  signal  and  the  trigger  would  result  in  phase  jumps  of  up  to  a  maximum  value 
of  7i  as  shown  in  Figure  4(a).  These  jumps  could  be  removed  by  dynamically  triggering  the  ADC 
using  an  electrically  tunable  TTF  signal  generated  from  a  narrow  band  (0.1  nm)  Fiber  Bragg 


Grating  (FBG)  as  shown  in  Figure  3.  A  reflected  optical  pulse  is  generated  whenever  the  source 
sweeps  the  FBG  reflection  wavelength.  An  optical  pulse  is  detected  by  the  detector  every  time 
the  laser  swept  1315  nm  (the  reflection  wavelength  of  FBG),  thus  the  frequency  of  the  FBG 
pulses  are  also  at  20  kHz  as  shown  in  Figure  4(a).  Figure  4(b)  depicts  the  oscilloscope  trace 
showing  the  TTL  trigger  signal  (2nd  trace),  the  signal  from  the  sample  (3rd  trace)  and  the  MZI- 
OC  signal  (4th  trace).  This  pulse  is  converted  into  a  TTL  signal  using  a  signal  generator  (Stanford 
Research  Systems).  The  TTL  signal  generator  generates  an  optical  pulse  which  can  be 
electrically  tunable  to  get  the  required  duty  cycle  as  shown  in  the  Figure  4(a).  By  triggering  the 
ADC  with  the  above  TTL  signal,  the  jitter  due  to  electronics  is  reduced  by  introducing  perfect 
synchronization  between  the  source  and  data  acquisition,  which  in  turn  reduces  the  phase 
variations:  standard  deviation  of  0.016  radians  for  512  A-line  scans  has  been  achieved  as  shown 
in  Figure  4(d). 

Regardless  of  the  above  removal  of  n  jumps,  the  MZI-OC  introduces  some  phase  noise  to  the 
system.  This  is  due  to  the  fact  that  zeros  and  peaks  keep  changing  with  every  laser  scan  due  to 
thermal  and  mechanical  vibrations.  Thus,  the  phase  measured  from  the  sampled  fringe  signal 
inhibits  high  phase  variations.  To  avoid  this,  the  MZI-OC  is  recorded  only  once  and  hence  the 
zeros  and  peaks  will  not  change  with  every  scan  and  the  recalibrated  signal  will  be  stable  in  the 
phase.  The  obtained  uncertainty  in  phase  is  very  low  and  is  experimentally  verified  in  the 
following  manner.  The  phase  shift  measured  in  water  as  a  function  of  time  for  different  glucose 
concentrations  and  corresponding  change  in  sample’s  refractive  index  as  a  function  of  glucose 
concentration  are  shown  in  Figure  5a  and  Figure  5b,  respectively.  Glucose-induced  changes  in 
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Figure  5:  (a)  Phase  shift  measured  in  a  clear  aqueous  solution  of  glucose  vs.  time  (or  scan  number);  and  (b)  change 
in  the  sample’s  refractive  index  as  a  function  of  glucose  concentration. 

the  refractive  index,  — ,  were  calculated  as  —  =  where  d  is  the  cell  thickness  and  —  is 

dC  dC  d  4jt  dC  dC 

cLti 

the  glucose-induced  changes  in  the  phase.  The  standard  deviation  for  each  —  data  point  was 

calculated  from  80  independent  measurements  of  the  same  glucose  concentration  (Figure  5b). 
Linear  fitting  of  the  data  points  using  a  linear  least  squares  algorithm  yielded  a  correlation 
coefficient  r>0.999  and  p<  0.01%.  The  phase  increased  in  equal  steps  indicating  that  the  phase 
measured  is  indeed  the  phase  introduced  by  the  change  in  concentration  of  the  glucose  solution 


with  low  uncertainty  (0.03  radians).  With  the  stabilized  phase,  the  refractive  index  change  as 
small  as  6.3  x  10"6  (through  1  mm  optical  path)  can  be  detected  using  the  PhS-SSOCT. 

Results  and  Discussions 


Initial  testing  in  clear  media 
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The  developed 
system  was  initially 
characterized  by  using  a 
cover  glass  slide  (0.16 
mm  thick)  placed  in  the 
sample  arm  of  the 
interferometer.  A  cover 
glass  slide  is  placed  at 
five  different  distances 
from  the  reference  arm 
and  imaged  without  and 
with  MZI-OC 

calibration  (shown  in 
Figure  6(a)  and  Figure 
6(c),  respectively).  Figure 
6(b)  and  Figure  6(d)  show 
the  corresponding  ID 
depth  profiles.  These 
data  clearly  demonstrate 
that  the  resolution  is 
degraded  without  proper 
calibration  and 

remapping  into  /c- space. 
Figure  6(c)  and  Figure 
6(d)  shows  that  after 
proper  recalibration  and 
remapping  into  k-space, 
there  was  no  degradation 
in  resolution  and 
maximum  imaging  depth 
of  6  mm  is  achieved. 


Figure  6:  2D  images  of  a  glass  slide  of  thickness  0.16mm  at  different  depths 
before  (a)  and  after  (b)  recalibration  and  remapping  to  equal  /c-space  and  (c),  (d) 
corresponding  1-D  depth  profiles.  A,B,C,D,E,F  are  the  A-line  scans  at  different 
depths. 


In  general,  the  path  difference  corresponding  to  the  3  dB  drop  in  the  1-D  profile  yields  the 
axial  resolution  of  the  system.  For  the  developed  PhS-SSOCT,  an  axial  resolution  of  10  pm  and 
a  transverse  resolution  of  12  pm  at  20  kHz  speed  of  an  A-line  scan  were  obtained.  Theoretical 

2fln 

axial  resolution  is  7  pm  (computed  as:  A z  =  — — — ).  The  mismatch  in  theoretical  and  measured 

values  is  attributed  mainly  to  dispersion  effects:  utilized  imaging  lens  (KBX025,  Newport 
Corporation)  disperses  the  ends  of  the  wavelength  scans  making  those  wavelengths  partially 
incoherent  which  effectively  reduces  the  wavelength  range  and,  thus,  increases  A z.  Transverse 
resolution  can  be  enhanced  by  using  a  high  numerical  aperture  lens  with  compromising  the 


imaging  depth.  However,  Y.  Verma  et.al.  recently  have  shown  that  the  use  of  a  tapered  single 
mode  fiber  could  increase  the  transverse  resolution  without  significant  decrease  of  the  imaging 
depth  [91], 

Performance  of  the  developed  system  was  further  evaluated  in  water  containing  gas 
microbubbles  of  different  diameters.  Water  was  injected  into  a  500  pm  flow-through  cuvette,  and 
bubbles  were  generated  by  introducing  different  pressures  using  a  peristaltic  pump  (Fisher 
Scientific).  The  beam  is  scanned  across  the  cuvette  as  shown  in  the  Figure  7(a).  The  amplitude  of 
the  interference  signals  in  the  time-delay  domain  was  recorded  from  the  cuvette.  Four 
characteristic  interferometric  peaks  were  observed  corresponding  to  the  interferences  between 
the  four  surfaces  of  the  cuvette  [an  example  is  shown  in  Figure  7(b)]. 

In  these  experiments,  the  optical  delay  is  calculated  as  a  function  of  dynamic  refractive  index 
(modified  by  presence/absence  of  microbubbles)  from  the  interferometric  peaks  that  are 
produced  by  the  reflection  from  the  inner  walls  of  the  cuvette  in  the  time-delay  domain  (between 
the  surfaces  3  and  2  or  B  in  Figure  7(b)).  The  phase  is  extracted  from  the  complex  Fourier 
transform  of  the  interference  fringes  and  monitored  at  the  interferometric  peak  that  corresponds 
to  the  self-interference  between  the  inner  glass  surfaces  3  and  2.  The  signal  from  surface  2  in 
Figure  7(a)  is  used  as  a  reference  arm  in  generating  the  self-interference  signal.  Phase-sensitive 
measurements  of  water  are  taken  before  and  after  injection  of  the  microbubbles.  Since  the 
bubbles  induced  changes  in  the  refractive  index,  the  diameter  of  the  bubble  from  the  phase  shift 
is  calculated  from  the  Equation  3: 


dn  =  ~  d<p,  (3) 

where  dn  is  the  refractive  index  change  introduced  by  the  bubble  (equal  to  0.33  for  an  air  bubble 
in  water),  /  is  the  diameter  of  the  bubble,  dcp  is  the  bubble-induced  phase  shift.  Phase  difference 
between  the  A  line  corresponding  to  the  center  of  the  bubble  (determined  as  a  maximal  possible 
change  in  the  phase)  and  the  A-line  outside  the  bubble  are  computed  to  quantify  the  diameters  of 
the  bubbles. 

Bubbles  are  generated  by  creating  a  pressure  difference  and  are  circulated  using  the  pump. 
Bubbles  of  diameters  greater  than  10  pm,  were  imaged  first  to  test  the  consistency  of  the  PhS- 
SSOCT  imaging.  The  images  were  taken  in  the  inter-interference  mode,  in  which  the  reference 
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arm  can  be  placed  in  such  a  way 
that  all  surfaces  of  the  cuvette 
(TYPE-40,  NSG  Precision  cells, 
Inc.)  can  be  shown  in  real 
distances  as  shown  in  Figure  7(a). 
Phase  stability  was  high  in  self- 
interference  mode  compared  to 
the  inter-interference  mode  as 
discussed  earlier.  Thus,  in  all  the 
results  presented  here,  the 
images  were  taken  in  inter¬ 
interference  mode,  and  the  phase 
response  was  measured  in  self¬ 
interference  mode. 
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Figure  7:  (a)  Orientation  of  the  cuvette  with  respect  to  the  beam  (b) 
corresponding  1-D  depth  profile  (A-  glass  thickness,  B-  optical  path 
thickness,  C-  A+B,  1,2, 3, 4  are  the  surfaces  of  glass  cuvette). 
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Figure  8(a)  shows  the  image 
of  a  500  pm  cuvette  fdled  with 
water  containing  a  bubble  of 
diameter  224  pm  (actual 
diameter).  A  high  display 
threshold  is  selected  so  as  to 
the  self-interference 
image.  The  1-D 
profile  of  the 
corresponding  self¬ 
interference  image 
and  the  temporal 
phase  response  are 
shown  in  Figure  8(b) 
and  Figure  8(c), 
respectively.  It 

should  be  noted  that 
the  imaging  and  the 
phase  response  were 
taken  one  after  the 
other  in  very  quick 
succession,  as  the 
imaging  was  done  in 
inter-interference 
mode  and  the  phase 
response  was  studied 
in  self-interference 
mode. 


Figure  8:  (a)  Image  of  a  cuvette  of  500  pm  thick  containing  a  bubble  of  224  pm 

(1,2, 3, 4  are  the  surfaces  corresponding  to  1,2, 3, 4  in  figure  4(a)  and  S  is  the  self-  The  reference 

interference  image)  (b)  corresponding  1-D  profile  showing  the  bubble  induced  shift  in  arm  js  placed  at  a 
peak  (c)  temporal  phase  response  at  665  pm  peak.  distance  SO  that  the 


whole  cuvette  can  be  imaged  without  any  negative  images  that  arise  due  to  the  symmetry  of  FFT. 
Therefore,  for  Figure  8(a),  all  four  bright  lines  correspond  to  all  four  surfaces  of  the  cuvette.  Each 
bright  line  in  the  Figure  8(a)  and  Figure  9(a),  correspond  to  each  surface  of  the  cuvette  (numbered 
1  to  4  in  the  Figure  7).  The  missing  part  in  the  surfaces  3  and  4  is  observed  only  in  the  case  of 
huge  bubbles  indicating  that  the  beam  has  refracted  so  that  the  reflection  from  the  surface  below 
the  bubble  is  not  coupled  back.  In  the  clear  aqueous  media,  the  maximal  phase  variations  were  as 
low  as  0.03  radians  after  a  5-point  averaging  (repeated  for  at  least  10  times),  which  implies  that 
any  microbubble  that  introduces  a  phase  shift  greater  than  0.03  radians  can  be  detected. 

Generally,  for  an  optical  path  difference  of  one  wavelength  the  phase  shift  would  be  2n  in  a 
homogenous  media.  Equation  3  describes  a  relationship  between  the  changes  in  the  refractive 
index  and  the  phase  of  PhS-SSOCT.  For  an  air  bubble  in  water,  the  change  in  the  refractive 
index  would  be  0.33  which  translates  to  a  minimum  bubble  size  of  2  pm  for  the  phase  to  be 
unwrapped  by  one  2n  jump.  Since  the  system’s  resolution  is  10  pm,  bubbles  with  diameters 
greater  than  10  pm  can  be  clearly  seen  in  the  structural  image,  as  shown  in  Figure  8(a).  When 
there  is  no  bubble,  the  optical  path  length  between  the  inner  surfaces  of  the  cuvette  is  665  pm 
(refractive  index  of  water  is  1.33,  so  500  pm  x  1.33  =  665  pm),  which  is  observed  as  a  peak  at 
665  pm  in  corresponding  1-D  depth  profile.  As  the  beam  interacts  with  the  bubble,  the  optical 
path  length  keeps  decreasing  until  it  reaches  the  center  of  the  bubble  and  increases  again  to  the 
original  value.  This  change  in  the  optical  path  length  is  reflected  as  a  shift  in  the  peak  in  the 
corresponding  1-D  profile  as  shown  in  the  Figure  8(b).  The  larger  the  bubble,  the  greater  is  the 
decrease  in  the  optical  path  length  and  hence  the  greater  is  the  shift.  The  number  of  2n  jumps  by 
which  the  phase  to  be  unwrapped  is  then  calculated  from  the  peak  shifted  in  the  1-D  profile.  This 
unwrapped  phase  is  then  added  to  the  PhS-SSOCT  phase  response  to  get  the  true  phase  response. 
For  the  clear  media,  each  depth  pixel  corresponds  to  a  one  way  physical  path  difference  of  5.88 
pm,  which  is  equal  to  2.94  “27i”  jumps.  Thus,  for  the  bubble  shown  in  Figure  8(a),  the  number  of 
2ti  jumps  would  be  69/2*2.94  =  113  as  the  peak  shifts  by  69  depth  pixels  [Figure  8(b)].  By 
plugging  this  true  phase  in  Equation  3,  the  size  of  the  bubble  obtained  is  219  pm.  The  actual 
diameter  of  the  bubble  measured  is  224  pm.  The  true  phase  and  the  PhS-SSOCT  phase  (phase 
before  adding  required  2n  jumps)  is  plotted  in  the  same  graph  with  two  different  scales  on  the 
right  and  left  of  the  Y-axis. 

Similarly,  several  bubbles  with  different  diameters  (52  pm,  94  pm,  160  pm  etc.,)  were  taken 
and  quantified.  The  obtained  error  ranged  from  0.19  pm  to  10  pm  which  can  be  attributed  to  the 
2n  ambiguity.  The  change  in  the  path  differences  less  than  10  pm  is  not  reflected  in  the  1-D 
profile  (due  to  limited  imaging  resolution  of  10  pm).  In  this  case,  the  number  of  27i  jumps  by 
which  the  phase  should  be  unwrapped  cannot  be  determined  and  can  take  any  integer  value 
between  1  and  5.  This  means  that  microbubbles  with  diameters  with  multiples  of  2  pm  (up  to  10 
pm)  would  show  the  same  phase  response.  For  instance,  if  the  number  of  27i  jumps  calculated 
from  the  phase  shift  is  60,  the  actual  number  of  2n  jumps  can  be  any  number  between  60  and  65. 
Likewise,  if  the  measured  diameter  of  a  bubble  is  120  pm,  then  the  actual  bubble  diameter  could 
be  120,  122,  124,  126  or  128  pm.  However,  this  ambiguity  can  be  resolved  by  implementing  a 
fast  real-time  continuous  unwrapping  algorithm  which  will  be  developed  in  our  future  studies. 
This  algorithm  would  acquire  multiple  phase  recordings  between  consecutive  changes  in  the  path 
difference  of  less  than  2  pm  and  continuously  unwrap  the  phase  information.  It  would  also  allow 
the  exact  quantification  of  microbubbles  that  are  beyond  the  imaging  capabilities  of  the  system. 


Figure  9(a)  show  an  example  of  small  bubbles  which  cannot  be  resolved  from  SSOCT 
structural  imaging.  These  bubbles  can,  however,  be  detected  by  the  PhS-SSOCT.  The  sizes  of 
the  bubbles  were  estimated  to  be  1.9  pm  (but  could  be  3.9  pm,  5.9  pm,  7.9  pm  or  9.9  pm  due  to 
the  ambiguity  discussed  above).  Therefore,  PhS-SSOCT  is  an  effective  device  for  ultra-sensitive 
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Figure  9:  Image  of  a  cuvette  of  500  pm  thick  with  (a)  very  small  bubbles  and  corresponding  temporal  phase 
response  at  665  pm  peak  (b)  of  (a). 

quantification  of  microbubbles  with  diameters  significantly  less  than  imaging  capabilities  of  the 
employed  system.  However,  the  error  due  to  ambiguity  can  still  be  decreased  by  increasing  the 
number  of  A-line  scans  with  in  the  same  frame  and  by  reducing  the  focused  beam  diameter.  For 
instance,  as  the  beam  propagates  across  the  bubble,  it  interacts  more  and  more  of  the  air/gas 
volumes  in  the  bubble  indicating  the  increase  in  diameter.  If  a  1 0  pm  bubble  is  scanned  in  5  steps 
of  2  pm  (the  difference  between  two  consecutive  A-line  diameters)  each  for  half  the  bubble,  then 
the  phase  would  be  needed  to  unwrap  only  once.  Hence,  no  2n  ambiguities  would  be  observed. 
To  achieve  that,  a  3  mm  transverse  length  should  be  scanned  with  1500  A-lines  as  opposed  to 
512  A-lines.  Thus,  to  obtain  1500  A-lines  in  the  same  time,  the  speed  of  the  system  should  be 
enhanced  by  about  3  times,  which  can  be  achieved  by  using  Fourier  domain  mode  locked  laser  as 
described  in  [92], 

Experiments  with  scattering  tissue-simulating  media 

The  detection,  imaging  and  quantification  of  microbubbles  in  scattering  media  were  also 
performed  on  a  similar  note.  1.54%  polystyrene  spheres  (PS)  were  used  to  simulate  the  scattering 
media  with  a  scattering  coefficient  of  100  cm"1  for  the  1324  nm  wavelength.  Figure  10(a)  depicts 
the  image  of  the  cuvette  when  filled  with  scattering  media.  A  high  dynamic  range  is  chosen  to 
identify  if  any  self-interference  lines  are  observed.  The  phase  is  monitored  at  the  peak 
corresponding  to  705  pm  path  length.  The  phase  variations  obtained  were  0.04  radians  (the  same 
phase  sensitivity  is  obtained  even  if  milk  is  used  as  scattering  media),  indicating  the  minimum 
size  of  the  bubble  that  could  be  detected  is  0.01  pm.  Figure  10(b)  shows  the  image  of  a  bubble 


with  an  actual  diameter  of  142  pm,  and  the  diameter  obtained  from  PhS-SSOCT  measurements 
was  143  pm  from  Figure  10(d).  Note  that  the  phase  response  in  Figure  10(d)  is  shown  for  few  A- 
lines  so  that  the  phase  response  could  be  clearly  displayed. 
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Figure  10:  Image  of  the  cuvette  with  scattering  media  (a)  without  any  bubble  (c)  with  a  bubble  and  their 
corresponding  temporal  phase  response  at  the  self-interference  peak  in  (b)  and  (c). 

The  bubbles  that  are  beyond  the  imaging  capabilities  of  SSOCT  are  shown  in  Figure  11.  The 
portion  circled  in  blue  in  the  image  is  the  place  where  the  phase  indicates  that  there  could  be 
three  bubbles  of  diameters  all  around  0.8  pm  (but  that  could  be  2.8  pm,  4.8  pm,  6.8  pm  or  8.8 
pm).  As  the  sizes  of  the  bubbles  are  smaller  than  the  focused  beam  spot  (~  25  pm),  they  act  as  a 
hindrance  to  the  beam  passage  and  create  a  shadow  on  the  other  interface  as  shown  in  Figure 
11(a). 

The  PhS-SSOCT  is  also  tested  for  detection  of  fast  moving  large  bubbles.  Figure  12  depicts 
the  M-mode  image  of  air  bubbles  and  air  gaps  that  encountered  the  beam.  The  sizes  of  the 
bubbles  can  be  predicted  from  the  depth  scan  (Y-axis)  and  the  speed  of  the  bubbles  from  the 
transverse  scan  (X-axis).  The  large  air  gaps  have  a  bright  line  at  the  end  of  the  cuvette  which  was 
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Figure  11: 

response. 


(a)  Cuvette  with  very  small  bubbles  and  (b)  corresponding  temporal  phase 


less  deep  than  when 
there  was  no  air  gap, 
indicating  the 

presence  of  some 
medium  that  has  a 
lower  refractive 
index  than  the 
scattering  media. 

The  moving 
bubbles  could 

introduce  a  shift  in 
phase  due  to  the 
Doppler  effect.  The 
maximum  phase  shift 
that  can  be  observed 
in  a  capillary  where 
the  maximum  blood 
flow  could  be  1 


mm/s,  is  0.2  radians  (computed  using  the  equation  3  with  p  =  85  and  r  =  0.15  ms). 


v  = 


A(J> 


2n<k>TCOs(/3) 


(4) 


where  V  is  the 
velocity  of  the  bubble 
obtained  from  the 
Doppler  image  of  the 
blood  flow,  k  is  the 
wave  number,  r  is  the 
time  taken  between 
two  A-lines  and  P  is 
the  angle  made  by  the 
velocity  vector  with 
the  laser  beam.  The 
error  in  quantifying 
the  micro-bubbles  that 

could  be  introduced  by  this  phase  shift  from  the  Doppler  signal  could  be  +/-  0.06  pm.  However, 
assuming  that  the  bubbles  and  the  blood  flow  at  the  same  velocities,  this  error  can  be 
compensated  by  measuring  the  phase  shift  due  to  the  Doppler  effect.  Our  future  studies  include 
the  Doppler  imaging  for  the  quantification  of  moving  micro-bubbles  in  a  blood  flow.  Further 
experiments  would  be  conducted  with  moving  bubbles  to  validate  the  compensating  factor. 


Figure  12:  Image  of  500  (im  cuvette  containing  fast  moving  bubbles. 


In  Vivo  Experiments 

Encouraging  results  were  obtained  from  preliminary  in  vivo  studies  with  mice. 
Microbubbles  used  for  these  in  vivo  studies  were  50%  Definity  microbubbles  in  saline  and  air 
microbubbles  generated  by  inducing  variable  pressures  in  the  syringe  before  injections. 


Continuous  imaging  of  the  mouse  tail  vein  at  a  location  2  to  3  cm  upstream  away  from  the  point 
of  injection  was  performed  before  injection,  during  the  injection  and  after  injection  of  both 
Definity  and  air  microbubbles. 


(a)  (b)  Before  animal  experiments, 

we  have  upgraded  the  laser  used  in 
the  system  and  the  re-developed 
sensor  shows  now  an  axial 
resolution  of  8  pm  and  transverse 
resolution  of  ~  15  pm,  A-line 
speed  of  30  kHz.  The  3D  imaging 
capabilities  of  the  system  are 
illustrated  in  the  Figure  13(a)  and 
Figure  13(b).  Figure  13(a)  shows  an 
example  of  a  2D  image  of  the 
blood  vessel  in  mouse  tail  and 
Figure  13(b)  shows  the  3D 
construction  of  the  vessel  from  its 
corresponding  Doppler  signal.  The 
2-D  images  required  to  reconstruct  the  3D  were  acquired  at  a  real  time  speed  of  40  frames  per 
second  with  the  frame  size  of  1600X500. 


Figure  13:  (a)  2D  SSOCT  image  of  a  mouse  tail  in  the  blood  vessel 
area  and  (b)  corresponding  3D  reconstruction  of  blood  flow  from 
Doppler  signal. 


(a) 


(b) 


Figure  14:  (a)  SSOCT  image  of  a  mouse  tail  in  the  blood  vessel  area 
and  (b)  corresponding  Doppler  signal  showing  blood  flow. 


100  pm  were  spotted. 


This  system  was  then 
utilized  for  detecting  Definity 
microbubbles  injected  into  blood 
vessels  of  mice.  Mouse  tail  vein 
was  detected  in  both  structural  and 
the  Doppler  image  as  shown  in  the 
Figure  14.  The  vein  is  then 
continuously  imaged  during  the 
injection  and  after  injection  of  the 
Definity  bubbles  for  12  minutes. 
Figure  15  shows  the  image  of  the 
vessel  9  minutes  after  injection  of 
Definity.  A  few  clusters  of  the 
Definity  bubbles  of  sizes  around 


Similarly  continuous  imaging  was  conducted  using  air  bubbles.  Figure  16(a)  shows  an 
example  of  air  microbubbles  of  various  sizes  (80  pm,  200  pm  etc.,)  detected  by  PhS-SSOCT  that 
were  induced  in  to  the  blood  stream  of  the  mouse  tail  vein.  Similarly,  a  bubble  of  size  56  pm  was 
spotted  in  the  tail  vein  is  shown  in  the  Figure  16(b).  Continuous  imaging  clearly  depicts  the 
movement  of  microbubbles.  Hence,  the  developed  sensor  was  capable  of  real  time  monitoring 
and  quantifying  of  microbubbles  in  animals  in  vivo. 


Conclusions 


Figure  15:  Cluster  of  Deflnity  microbubbles 
spotted  in  real-time  while  circulating  in  the 
tail  vein. 


Under  the  sponsorship  of  ONR  YIP,  we  have 
initiated  very  exciting  research  on  development  of 
novel  class  of  biosensors  for  detection  of 
decompression-associated  illnesses.  We  demonstrated 
that  the  developed  PhS-SSOCT  is  capable  of  imaging, 
detecting  and  quantifying  of  large  and  small 
microbubbles  in  clear  and  scattering  media  as  well  as  in 
blood  microvessels  in  vivo.  The  images  shown  were 
taken  using  the  conventional  SS-OCT  and  the  phase 
results  were  obtained  using  common  path  SSOCT.  The 
results  suggest  that  small  microbubbles  with  diameter 
beyond  imaging  capabilities  of  the  system  can  be 
detected  and  quantified  using  the  PhS-SSOCT. 
Potentially,  microbubbles  of  diameters  as  small  as  0.01 
pm  (that  introduce  a  phase  shift  of  0.03  radians)  can  be 
detected  and  quantified  with  this  method.  Our  future 
studies  will  focus  on  the  development  of  the  effective 
phase  unwrapping  algorithm  that  will  quantify  the 
microbubbles  without  any  2ji  ambiguities.  To  our 
knowledge,  OCT  was  utilized  for  the  first  time  to  detect 
microbubbles  in  animals  in  vivo.  Moreover, 
microbubbles  of  sizes  ~  50  pm  were  detected  for  the 
first  time  in  animals  in  vivo.  Thus,  from  these  results  it 
can  be  inferred  that  the  SSOCT  could  be  built  into  a 
real  time  useful  diagnostic  tool. 

Publications  resulted  from  ONR  grant  support 

The  described  results  were  published  in  2  original 
and  10  contributed  peer-reviewed  manuscripts  [93-104] 
and  one  book  chapter  [105]. 
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